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1.  INTRODUCTION 


The  results  of  past  studies  on  Lounded-wave  simulator  field  environ¬ 
ments  in  the  working  volume  (Refs.  1  through  9)  are  strictly  applicable  to 
ALECS  and  ARES  which  do  not  have  a  wood  stand.  The  wood  stand  (Fig.  1) 
is  an  integral  part  of  the  ATLAS  L  (Trestle)  and  its  presence  will  not  onl; 
affect  the  field  in  the  working  volume,  hut  also  modify  the  responses  of  trie 
test  object.  The  former  effect  is  deemed  more  -significant  than  the  latter. 
Numerical  studies  of  the  for"’  r  ,  have  been  accomplished  in  Reference  10. 

However,  analytical  studies  are  ex,  ruiuclv  invaluable  to  an  overall  understanding 
of  the  ATLAS  I  field  environment  in  the  working  volume,  and  have  never  been 
conducted. 

The  objective  of  this  report  is  to  develop  a  simple  electromagnetic 
model  for  the  ATLAS  I  wood  stand  that  will  yield  a  simple  engineering  analytic 
form  to  describe  accurately  the  ATLAS  I  field. 

In  Section  II,  formulas  will  be  developed  for  the  field  distributions 
and  dispersion  relations  of  the  TE  (it  .nsverse  Electric)  surface-wave  modes 
supported  by  dielectric  i  -  .  Section  III  is  devoted  to  deriving  numerical 
values  appropriate  for  the  ATLAS  I  simulator  based  on  the  formulas  in  Section  II. 
In  Section  IV  the  numerical  values  wi I  !  be  compared  with  the  ATLAS  I  field  mapping 
data  (Ref.  11)  and  simple  engineering  analytical  formulas  for  the  ATLAS  I  simula¬ 
tor  fields  will  be  obtained.  Finallv,  in  Section  V  a  summary  of  important  results 
will  be  given. 


5 


DISPERSION  RELATION 


1 


Given  a  three- lave  red  medium  (Fig.  2),  the  problem  at  hand  is  to  find 
the  propagation  constant  of  the  TE  surface  wave  (H^,  ,  E  ).  The  direction 

of  propagation  is  along  the  z  axis  and  there  is  no  variation  of  the  field  in 
the  y  direction.  The  time  variation  of  exp(jwt)  is  assumed  and  suppressed. 

To  find  the  propagation  constant  one  starts  with  Maxwell's  equations 


V  x  E  =  -  ju;uH 


7  x  H  = 
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Equation  6  will  be  solved  for  each  region  shown  in  Figure  2. 

1.  Region  0  (t  =  e^,  x  >_  a/2) 

Assume  there  is  no  propagation  in  the  x  direction  (decaying  wave). 
The  solution  of  Equation  6  for  E^^  can  be  written  in  the  following  form: 
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Equation  6  cakes  the  form 
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3.  Region  2  t  e 
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Equation  <>  is  of  tin-  form 
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the  solution  of  which  with  a>>  pro]  mi;  on  in  tin  -  direction  is  given  by 
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Matching  the  wave  impedances  across  the  interfaces  one  gets 
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Substituting  the  field  components  given  in  Equations  8,  12,  14,  17,  i9  and 
into  Equations  23  and  24  one  obtains  the  dispersion  relation 


,2 

1  +  a  A 

tanh"  ( A.  a/2)  +  2  — tt-~ — r  tanh(A  a/2)  +  1  =  0 

1  a  +  A  )  i 


(25) 
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which  can  be  written  in  the  following  form 


tanh(A^a) 


v1(',+  v  ) 
12  o 


A ,  +  A  A0 
1  o  2 


(26) 


tanh 


(..A  ) 


(27) 


Equation  27  is  the  final  form  for  the  dispersion  relation  of  the  problem, 
and  will  be  solved  for  (=  B/k  )  in  the  next  section. 


One  could  have  del  ined  -  k  /  -  in  Equation  15  instead  r-’ 

n  1  °  lr 

'  ~  £^r-  Then  the  ii.-perbolic  functions  sinh,  cosh  and  tank: 
in  Equations  14,  17,  25  and  do  could  have  become  trigonometric  functions 
sin,  cos  and  tan  after  making  the  following  replacements 


1  i  ‘  1 

sinh  (AjX)  -  j  sin  (  •  x) 
cosh  (  1 x )  •  cos  i  jx) 
tanli  ( 1  !  x)  •  j  tar.  ( A  ^  x) 
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III.  PROPAGATION  CONSTANT  AND  WAVK  1 MPKDANCK 


To  find  the  normalized  propagation  constant  £  for  given  values  of 

t',  ,  e„  and  k  a  the  dispersion  haunt  ion  27  has  to  he  solved  numerical!'.', 
lr  2r  o 

Before  presenting  the  numerical  solution,  the  ranges  of  solution  for  real  l 

under  the  condition  of  e,  >  e,  >1  are  analvzed. 

1  r  _  r  — 

In  the  range  whore  l.~  1 ,  t  he  left-hand  side  of  Equation  27  is  purel  / 

imaginary,  but  the  right-hand  side  of  this  equation  is  real.  This  means 

t 

that  no  roots  are  possible  for  i,“~  <■  1. 

In  the  range  where  >  ’  i ,  '  !>  •  left-hand  side  of  Equation  27  is 

still  purely  imaginary  but  the  right-hand  side  of  this  equation  is  complex. 

In  order  for  this  equation  to  have  roots, the  real  part  of  the  right-hand 
side  must  be  equal  to  zero.  A  little  algebra  shows  that  this  Is  not  possible. 
Thus,  no  roots  exist  within  this  range. 

2 

In  the  range  where  >  5  i  e2r’  koth  the  Itft-  and  right-hand  sides 
of  Equation  27  are  purely  imaginary.  This  means  that  real  roots  for  £  may 
exist  in  this  range. 

2 

In  the  range  where  £  2.  cir’  !  ho  right -hand  side  of  Equation  27  is 

negative  real  and  left-hand  side  is  positive  real.  Hence,  no  roots  are 
possible  in  this  range. 

The  above  simple  analysis  shows  that  the  real  normalized  propagation 
constant  £  is  limited  to  the  range  >  £  >  /ejr .  In  the  following,  two 

cases  will  be  considered,  namely  e0r  =  1.04  and  =  !• 

1.  c,  =  1.04 
2r 

From  Reference  12,  the  dimensions  and  spacings  of  the  wooden  struts, 

the  effective  dielectric  constant  ol  the  region  below  the  wood  platform  is 

estimated  to  be  about  1.04.  Figure  3  shows  £  of  the  first  TE  surface-wave  mode 

versus  k  a  for  e,  =1.04  and  e,  =  4,6,10. 
o  2r  lr 

An  Inspection  of  Figure  3  and  Equation  27  reveals  that  no  real  £,  i.e., 

no  propagation,  is  possible  until  some  critical  value  of  k  a  is  reached.  Let 

o 

this  frequency  be  called  the  cutoff  frequency  of  the  surface  wave. 
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Obviously,  the  cutoff  frequency  can  be  obtained  from  Equation  27  with 
2 

£  =  e2r>  that  is. 


tanh  (k^n)a  )  = 


/cl -  c,  /cT  -  i 
_ 2r  lr  2r 

e,  -  c., 
lr  2r 


which  gives 


.  (n  -  l)n  +  arctan  (  /(e..  -  1)/  U,  -  c,  )  j  o 

fln)  _ _ _ _ 1. - ±£ - — - 61 -  x  3  x  io° 

c  2va/c.  -  c0 

lr  2r 


Here,  is  the  cutofi  frequency  of  tn-.  mode  of  the  surface  wave.  The 

c 

first  three  cutoff  frequencies  are  given  in  Table  1. 

TABLE  1.  FIRST  THREE  CUTOFF  FREQUENCIES 


n 

£, 

lr 

1 

2r 

— 

a  (mi 

f(n)Oniz) 

C 

- 

k(n)a 

c 

i 

4 

1.04 

1 

0.0673 

2 

4 

A  .  W  t 

1 

o.  • 

1.8933 

3 

4 

1.04 

1 

177.6 

3.7193 

(1). 


Figure  4  shows  the  normalized  cutoff  frequency  (k^  a)  versus  c2r  f°r  Eir~  ^ 


2.  e2r=l 


When  e„  =1,  Equation  27  becomes,  for  modes  with  H  (-x)  =  H  (x), 
zr  x  x 


/2)‘  - 


tanh  \  kQa, lr 


and  becomes,  for  modes  with  H^C-x)  =  -  H^fx) 


^2-  1 
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(28) 


lr 


)th  (k  a  / f  2  -  ■ .  H  )  =  -  — 

\  o  lr  /  r 2 - 


lr 


-10. 
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Figure 


The  normalized  propagation  constant  £  has  real  values  under  the  condition 

e  >  1  at  any  frequency.  Figure  5  shows  f.  of  the  first  TE  surface-wave  mode 

versus  kQa  for  various  e^.  It  clearly  shows  that  there  is  no  cutoff  frequency 

for  e-  = 1. 

2r 

In  Figure  5,  the  C-value  versus  kQa  for  e2r  =  1.04  is  also  given.  It 

is  observed  that  the  difference  in  ^-values  for  -  1  and  t2r  =  l*04  is 

generally  negligible  except  at  very  small  k^a.  Thus,  only  the  simpler 

results  with  e„  = 1  will  be  used  in  the  following  discussion  and  for  the 
2r  * 

comparison  with  the  field  mapping  data. 

After  obtaining  the  ’  "'ncdances  ^  and  ^  can  be 

calculated  from  the  following  two  equations: 


k  -  1 
o 


42- 1 


6  ~  "  °/ko  C 


which  are  plotted  in  Figures  6  and  7  for  =  I  anc*  various  e^. 

In  the  next  section,  the  above  results  will  be  used  to  compare  with 
the  ATI  AS  I  field  mapping  data. 


Since  Table  1  shows  that  the  first  cutoff  frequency  is  about  3.2  MHz  for 
4  and  t2r  =  1.04,  one  may  question  the  validity  of  neglecting  the 
lower  medium  with  E2r  =  1.04.  However,  at  this  frequency  the  effect  of 
the  ground  is  no  longer  negligible.  When  the  ground  is  properly  taken  into 
account,  there  will  be  no  cutoff  phenomenon. 
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Normalized  of  the  first  TE  surface-wave  mode 

1  y  z  ' 

versus  k  a  for  e„  =1  and  e,  =4,6,8,10. 
o  2r  lr  ’  ’  ’ 


i  i 


Figure  7. 
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IV.  ANALYTICAL  REPRESENTATION  OF  ATI  AS  I  FIELDS 


Hie  first  TE  surface-wave  mode  described  in  the  previous  sections 
will  be  used  to  represent  the  field  distributions  within  the  TRESTLE 
simulator. 

In  Section  III,  curves  for  the  ratio  of  the  electric  to  magnetic  field 

(i.e.,  the  wave  impedances)  of  the  first  TE  surface  wave  versus  frequency 

were  obtained.  The  impedance  curves  (without  normalization)  for  =  4, 

10,  and  a  =  0.5,  1,  1.5m  are  plotted  in  Figures  8  and  9  in  the  log-log 

scale.  In  the  figures,  the  imr  'dinc'-s  deduced  from  the  ATLAS  1  field 

mapping  data  (Ref.  1 1)  are  also  given.  Good  agreements  are  observed 

between  the  field  mapping  data  and  the  analytical  results  based  on  the 

first  TE  surface-wave  mode  when  e ,  -  4  and  a =  1 .  (Also,  see  the  final 

lr 

remark  on  pages  29  and  30) . 

In  Figures  8  and  9,  all  the  available  field  mapping  data  have  been 
used  to  obtain  the  impedances  except  for  test  points  17  and  21  whose 
impedances  are  expected  to  be  approximately  equal  to  those  of  test  points 
13  and  22  (Fig.  10).  Figure  10  shows  all  the  test  points  that  are  in 
Reference  11. 

The  good  agreement  shown  in  Figures  8  and  9  gives  one  the  confidence 
in  using  the  first  TE  surface-wave  mode  for  desci ibing  the  ATLAS  I  fields. 
The  next  step  is  the  determination  of  the  constant  A  in  Equation  8.  To 
this  end,  the  field  mapping  data  of  test  point  2  are  used.  One  typical 
set  of  the  frequency-domain  curves  of  and  at  this  test  point 

is  given  in  Figure  11.  The  frequency  dependences  of  and  are 

almost  the  same,  as  they  should  be  according  to  Equation  11,  for  the 
frequency  range  where  kQa  <  0.2  (frequency  <  10  MHz).  The  asymptotes 
are  drawn  in  Figure  11  in  broken  lines  leading  to  the  following  form  for 
the  constant  A: 

Z  H 

A  (1+  st^Kl  +  st^)  (31) 
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Figure  11. 


Frequency-domain  curves  of  measured  and  at 
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test  point  2  and  their  asymptotes  (Ref.  11). 
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whore 


H  =  6  »  10  ^A/ia-Hz 
o 

V3*  10“7S 

t2  =  8  x io"9s 


In  obtaining  Equation  31,  one  has  assumed  >^x  =  0  at  test  point  2,  and 
ignored  the  notch  effect. 


From  Equations  8,  11  and  31, 
test  point  2  are 


th  •  time-domain  E^(t)  and  H^(t) 
y  x 


at 


(t)/Z 

o 


(32) 


which  is  plotted  in  Figure  12.  The  agreement  with  field  mapping  data  at 

late  times  is  excellent.  But  there  is  considerable  difference  at  early 

times.  The  difference  is  attributable  to  the  notch  existing  in  the  pulser 

voltages.  To  account  for  the  notch  effect,  one  may  subtract  a  term  from 

Equation  32.  From  Figure  11  one  can  see  that  in  the  frequency  domain  the 

term  to  be  subtracted  should  behave  as  [l/s  I  when  |sj  is  large,  and  should 

have  a  double  pole  at  s  3  MHz.  Such  a  terra  has  a  time  variation  of  the 

form  t  exp(-t/t  )  with  t  =6x10  s.  Thus, 
o  o 


H<°> 

X 


(t)/z0= 


1.8t 

t 

o 


(33) 


where  the  coefficient  1.8/t^  is  chosen  in  such  a  way  that  the  best  agreement 
between  Equation  33  and  the  field  mapping  curve  can  be  obtained  (except  for 
the  prepulse  region.  Figure  12). 

From  Equations  12  and  33  the  z-component  of  the  ATLAS  I  magnetic  field 

can  also  be  estimated.  Generally,  it  will  involve  solving  Equations  28  and  9 

to  obtain  as  a  function  of  id  and, subsequently, inverting  complicated 

Fourier  (or  Laplace)  integral.  However,  if  one  is  only  interested  in  the 

late-time  behavior  of  H^(t)  where  the  high-frequency  part  of  the  spectrum 

Z  (0) 

is  not  important,  a  simple  expression  for  (t)  can  be  obtained  in  the 
following  manner:  From  Equation  12,  one  has 
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Figure  12,  Time-domain  curves  of  and  E ^ /7.  at  test  point  2  from  field 

x  y  o  r 

mapping  data  (Ref. 11)  and  analytical  representations  (Equations 
32  and  33). 
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(34) 


„(0),  ,  ,  Ao  '  (0) ,  .  _  .  Ao  '(0),  v 

Hz  (u)  - " J  EZ  (<j)  - j  r  Hx  (u) 

0 


u(c  -  1)«  .(0) 

j  — ±j -  Hx  <„>  »„.«!> 


Thus, 


z  2c  3t  x 


(elr-l)a  H 


2c 


H  r  i 

(tl  *t2) 


-t/t1  -t/c2 

e  -  —  e  + 

C2 


H(l-r)e  t/to]  <35) 


(for  t  >  t^) 


Equation  35  is  plotted  in  Figure  13  where  a  typical  field  mapping  curve 
is  superimposed.  The  estimated  H^°\t)  resembles  the  field  mapping  curve, 
although  relatively  low  in  magnitude.  The  under-estimate  of  the  late-time 
H^^(t)-value  is  probably  due  to  the  following  reasons: 

1.  The  sensor  used  in  the  field  mapping  test  did  not  have  an  accurate 
response  at  low  frequencies,  or,  more  specifically,  gave  an  over¬ 
estimate  at  the  lew  frequency  region  (Fig.  9). 

2.  The  fact  that  the  wood  platform  is  of  finite  extent  is  not  taken 
into  account  in  the  theory. 

It  should  be  noted  that  Equations  32,  33,  and  35  are  derived  for 

a  field  point  not  too  high  above  the  wooden  platform.  Thus,  they  are  valid 

only  at  field  points  where  x  <  1  for  the  important  spectrum  range. 

From  Figure  7,  it  is  observed  that  X  a  =  0.06  for  k  a  =  0.2  and 

o  o 

e^r  =  4.  This  means  that  for  a  =  lm  and  frequency  =  10  MHz,  the  decaying 
distance  D,  defined  by  XqD  =  1  is  approximately  given  by 

D  -  16.5  m  (36) 

Also,  0  is  proportional  to  (frequency)-2  at  lower  frequencies.  It  is 
therefore  reasonable  to  conclude  that  Equations  32,  33  and  35  are  satis¬ 
factory  representations  for  the  ATLAS  I  fields  up  to  as  high  as  15  meters 
above  the  wooden  platform. 


27 


0  500  1000  1500  2000 

Uns) 


Figure  13.  Time-domain  curves  of  at  test  point  2  (Ref. 11)  from  field 

mapping  data  and  analytical  representation  (Equation  35) . 


Although  =  4  and  a  =  lm  have  been  selected  for  the  above  analysis, 
other  values  that  show  good  agreements  in  Figures  8  and  "  can  be  used  as 
well  (e.g.,  e^r  =  6,  a  =  0.5m).  The  reason  for  selecting  Elr  =  4  and  a  «  lm 
is  that  they  are  closer  to  the  actual  situation.  However,  no  matter  what 
values  are  used.  Equations  32,  33  and  35  still  hold  true  while  the  D-values 
(i.e.,  Equation  36)  will  vary  somewhat. 
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V.  SUMMARY 

The  following  simple  analytical  expressions  have  been  obtained  for 
the  fields  above  the  wood  platform  of  the  ATLAS  1  simulator: 


E(0)(t)  =  7.5  x  104  (  e~3. 3*10  ”t  _  e-1.2xl08t  _  3  x  ^  e-l  W 

y 


BW(t)  =  2.5  x  10 
x 


~4(e-3- 


3xl06t  .  e-l-2  L08t  _  3  x  107t  e-1.6xl0?t )  2 


B(0)  (t)  =  4.2  x  10~6  (  e-3-3xl°  t-36e  1,2x10  1  +  9(1  -  l.6xlo?t)e 
z 


-1.6xl07t] 


.  9 
Wb/m 


where  B^(t)  is  valid  only  for  t  >  3  x  10  ?s  .  These  fields  can  be  compared 
with  the  criteria  EMP  fields  given  by  (Ref.  13) 


E(t)  =  5.24  x lu’  \  e 


4 (  4  x  l06t  -5  x l08t  ) 

\  e  -  e  / 


B(t)  =  1.75  x  10~4  (  e-4  X  10  -  e-5  x 10  c  )  wb/m2 

The  term  that  corresponds  to  the  notch  has  a  double  peak  on  the  negative 
real  axis  of  the  s-plane.  This  double  pole  lies  between  the  two  single  poles 
that  correspond  to  the  double  exponentials. 
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